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ABSTRACT 

We present combined optical/near-IR photometry (BVIK) for six open clusters - M35, 
0\ ; M37, NGC 1817, NGC 2477, NGC 2420, and M67. The open clusters span an age range 

from 150 Myr to 4 Gyr and have metal abundances from [Fe/H] = -0.3 to +0.09 
dex. We have utilized these data to test the robustness of theoretical main sequences 
constructed by several groups as denoted by the following designations - Padova, 
• Baraffe, Y^, Geneva, and Siess. The comparisons of the models with the observations 

(3 ; have been performed in the [My, {B - V)o], [My, {V - I)o], and [My, {V - K)o] 

If^ . colour-magnitude diagrams as well as the distance-independent [{V — K)o, {B — V)o] 

■ and [{V — K)o, {V — I)o] two-colour diagrams. We conclude that none of the theoretical 

I models reproduce the observational data in a consistent manner over the magnitude 

CO . and colour range of the unevolved main sequence. In particular, there are significant 

' zeropoint and shape differences between the models and the observations. We speculate 

, that the crux of the problem lies in the precise mismatch between theoretical and 

^ i' observational colour-temperature relations. These results underscore the importance 

of pursuing the study of stellar structure and stellar modelling with even greater 
intensity. 



o 



^2 ' Key vifords: Hertzsprung-Russell (HR) Diagram - open clusters and associations: 



general - stars: evolution - stars: fundamental parameters - stars: late-type 



C3 1 INTRODUCTION in the HR Diagram (e.g. Baraffe et al. 1998; Bergbusch & 

VandenBerg 2001; Yi et al. 2001). Aside from the uncertain- 

Witn the possible exception or the white dwari cooling se- ^- ■ ^ j_ • j_i ■ j_ • j_ • n j_i 

, ties inherent m their trigonometric parallaxes, the primary 

quence, the mam sequence is perhaps the best understood , , , r ■ . i . . . . . i ■ i • . i 

r 11 1 ^ r 1 , drawback or using these stars to test the isochrones is the 

phase or stellar evolution. Conversion oi hydrogen to helium , , , ,i ■ , r ,i • i-rr- ■ . n- 

need to correct their colours tor their dirtermg metallicities. 

occurs m the stellar core representing a stable and long last- 

, r 1 i- A 1 1 ii • Another method that has been used to test the theo- 

mg phase ot evolution. As we proceed along the mam se- 

j_i li-r retical isochrones is exemplified by the work of Andersen et 

quence trom high mass stars to low mass, the production OI i\t j •• c i u /inn/i\ 



energy in the core transitions from the CNO bi-cycle to the 
proton-proton chain, while, at the same time, the primary 
energy transfer mechanism in the core goes from convection 
to radiative diffusion. In addition, the outer convection zone 



al. (1991), Nordstrom & Johansen (1994), and most recently 
by Lastennet & Valls-Gabaud (2002, and references therein). 
These authors have determined the masses, effective temper- 
atures, and luminosities of individual stars in nearby binary 



, • 1 i-i j_i j_ ' 1 r 11 j_- systems and compared their HR diagram locations to the 

deepens progressively until the star becomes tully convective j f- a 

11 r,ora^ /n xr i 1 i f^r,o\ A 11 I- xi r i predictious of theoretical models. While the method has a 

below ~O.35M0 (Baraffe et al. 1998). All ot these features , r 

. 1 „ . 1 ,1 number oi intrinsic advantages, its primary shortcoming is 

are modelled with varying degrees ot success m modern stel- 

1 . , 1 1 , . 1 the obvious requirement that stars must be m binary sys- 

lar structure and evolution codes. 

mi 7 r M j_ j_ j_ 11 11- tems and in relative proximity to the Sun to be studied in 

ihe deqree or success we attribute to stellar models is , , j 

detail 

primarily based on how well they match the observations. 

Comparisons such as this have been performed in a number ^ ^^'""^ "^^^hod used to test theoretical stellar mod- 

of ways. One method that has been used in the past involves ^1^' ^""^ ^^e one we are utilizing in the present work, in- 
comparing isochrones to the locations of nearby field stars ^ol^^^ comparing the theoretical isochrones with duster 

colour-magnitude diagrams (CMDs) and two-colour dia- 
grams (TCDs) in various passbands (e.g. Castellani et al. 
* aaron@astro.ufi.odu 2002; von Hippel et al. 2002; von Hippel, Sarajedini, & Ruiz 

t ata@astro.ufl.edu 2003). Since all stars are of the same age and metal abun- 
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dance, this approach is less uncertain than comparing with 
individual field stars; however, the reddening is still required 
to perform the comparisons in the TCD (though the simi- 
larity in the slopes of the reddening vector and the model 
loci mitigates this somewhat), and both the distance and 
reddening are needed for the CMD comparisons. With the 
veritable explosion of recent high quality CCD photometry 
for open clusters, this technique is becoming increasingly 
more powerful as we will demonstrate in the present work. 

This paper is organized as follows. Section 2 presents 
a discussion of the optical and IR photometry. We briefly 
compare and contrast the various stellar evolution models in 
Section 3. Section 4 represents the primary scientiflc results 
of the present investigation; herein, we discuss the successes 
and failures of the models in terms of how well they match 
the two-colour diagrams and colour-magnitude diagrams of 
open clusters. Finally, we summarize our results in Section 
5. 



2 OBSERVATIONAL DATA 

The optical photometry for the six open clusters considered 
herein is taken from the following sources: M35, NGC 2420, 
NGC 1817 (WIYN Open Cluster Study, Deliyannis et al. in 
preparation), M67 (Montgomery, Marschall, & Janes 1993), 
M37 (Kalirai et al. 2001), and NGC 2477 (Kassis, et al. 
1997). To complement the optical data, we utilize JHKs 
infrared photometry available from the Second Incremen- 
tal Data Release of the 2MASS Point Source Catalog^ For 
all clusters in our sample we use the same criteria for the 
2MASS data retrieval. First, source brightness is limited to 
6th magnitude or fainter due to the possible saturation of 
bright sources (Carpenter 2000). Second, we have initially 
chosen all cluster samples to have a radius of 30 arcmin. 
The radius is then reduced to as small as 10 arcmin so as 
to isolate cluster members and minimize fleld star contam- 
ination. Lastly, the 2MASS catalogue provides a read flag 
(rdjlg) which indicates how photometry of each star was 
measured. We have excluded from our sample any sources 
marked with a read flag of zero in any one band since this 
imples that the star was not detected in that band and the 
magnitude reported is therefore only an upper limit. We 
note that a majority of source magnitudes come from point 
spread-function fitting (rd_flg — 2), however, we have chosen 
to include stars with aperture photometry (rdJlg = 1) so as 
to allow the inclusion of brighter cluster stars. 

The 2MASS program uses a K-short (Ks) filter for their 
photometry, which we have chosen to convert to the Bessell 
& Brett (1998, hereafter BB) K-band using the transfor- 
mation equations derived by Carpenter (2001) adapted as 
follows: 

(J-Kbb) = [(J-/s's)-(-0.011±0.005)]/(0.972±0.006)(l) 
and 

Kbb = [ifs~(-0.044±0.003)]-(0.000±0.005)(J-74'sB).(2) 
The BB photometric system has been adopted by a majority 

^ http://irsa.ipac.caltech.edu/ 



of the stellar evolutionary models analysed in the present 
paper. 

We combine optical and infrared photometry for stars 
in each cluster that are common to both data sets. All pho- 
tometry files have RA and Dec available for each star. To 
combine the data, we calculate the spatial offset between 
each star in the IR dataset with every star listed in the cor- 
responding optical dataset. When the spatial offset is mini- 
mized to within a given threshold, the two stars are consid- 
ered a match. The threshold distance is chosen based on the 
combined RA and Dec errors and turns out to be 2.8 arcsec 
for all clusters. 

Twarog et al. (1997) have compiled a database of infor- 
mation for 76 open clusters, for which they provide inter- 
nally consistent values of reddening, distance modulus, and 
metallicity. All of the clusters discussed in this study are in- 
cluded in the Twarog et al. (1997) catalogue, so we choose to 
adopt their values for the reddening and metallicity. Cluster 
ages are taken from Grocholski & Sarajedini (2002). In the 
case of M35, its age is calculated using the method described 
in Grocholski & Sarajedini (2002), which involves adopting 
ages from the WEBDA^ database and offsetting them to the 
age scale of Sarajedini (1999). We note that while we have 
tried to make the cluster ages as internally consistent as pos- 
sible, small variations in age have no effect on the location of 
theoretical isochrones along the unevolved MS. Table 1 lists 
these cluster parameters along with the filter passbands of 
the available photometry. 

Since we are primarily interested in comparing theoret- 
ical isochrones to cluster photometry in a variety of pass- 
bands, we have chosen to determine the distance moduli in 
the following manner. Given the adopted reddenings in Ta- 
ble 1, we shift each set of isochrones (described in Sec. 3) 
vertically in the {V — I)o CMD until the model main se- 
quence at Mv=6-0 fits that of the clusters. Since M37 does 
not have I-band data, we utilize {V — K)o for the MS fitting. 
We have avoided using {B — V)o for this purpose because it 
is more sensitive than {V — I)o to metallicity variations at 
temperatures above ~ 4000K (My ~ 8.5) which could result 
in erroneous distance determinations in cases where the clos- 
est available isochrone metallicity differs somewhat from the 
actual cluster abundance (Baraffe 2003, private communica- 
tion). The results, along with the distances given by Twarog 
et al. (1997), are given in Table 2. In this way, we can be 
sure that the models fit the main sequences in one colour at 
a given point allowing us to investigate how well they fit in 
other colours at other points along the MS. 

To correct for interstellar extinction, we adopt the red- 
dening law derived by Cardelli, Clayton, & Mathis (1989), 
which, using their value of Rv ~ 3.1 and taking into ac- 
count the central wavelength of each filter, we find the fol- 
lowing relations for the near-IR passbands: Aj = 0.282Ay, 
Ah = O.lSOAv, Ak = O.llGAv. For observations of M 37, 
Ab ~ lAlAv and for all other optical photometry, Ajj — 
1.57Av, Ab = 1.35Av, Ar = 0.839Av, and Ai = 0.560Av. 

In Figs. 1 and 2 we present [Mv, {B - V)o] and [My, 
{V — I)o] CMDs, respectively, for all clusters using the full 
optical data and in Fig. 3 we show [Mv, {V — K)o] CMDs 
using the combined optical and infrared datasets. The clus- 

^ |http : // obs www .unlge . ch / webda / webda. html | 
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ters are plotted in order of increasing age from upper left to 
lower right. We remind the reader that M 37 does not have 
I band observations available and is not included in Fig. 2. 



3 THEORETICAL MODELS 

We now introduce and briefly compare the five sets of the- 
oretical isochrones that are used in the present analysis - 

Padova (Girardi ct al. 2002)^, Baraffc (Baraffc ct al. 1998)", 
(Yi et al. 2001)^ Geneva (Lejeune & Schaerer 2001)'', 
and Siess (Siess, Dufour, & Forestini 2000) While we do 
not intend this to be a comprehensive review of the stellar 
models, we point out some of the similarities and differences 
between the input parameters of the five groups as listed in 
Table 3. 

With respect to the opacities, all of the groups use ver- 
sions of the OPAL opacities (e.g. Iglesias & Rogers 1996) for 
high temperatures supplemented by the Alexander & Fer- 
guson (1994, AF94) tables for low temperatures. The only 
caveat is that the Geneva models also consider the Kurucz 
(1991) low temperature opacities. We see that all five groups 
utilize different sources for their equations of state. Three of 
the groups use similar values for core overshoot in all of their 
models, while Siess has only included overshoot for their so- 
lar metallicity models. Baraffe has chosen to not include core 
overshoot because the focus of their work is low mass stars, 
which, for masses below ~ O.35M0, are fully convcctivo so 
an overshoot parameter is not necessary. Also, for low mass 
stars above this limit, a convective core is only present for a 
small fraction of time during their evolution, so inclusion of 
overshoot will not have a significant effect on these objects 
(Chabrier & Baraffe 1997). Mixing length parameters are 
similar for all groups, varying by less than 20%. It is impor- 
tant to keep in mind that the Baraffe models are available 
with three mixing length values, a = 1.0, 1.5, 1.9. Setting the 
mixing length equal to 1.9 is required to match observations 
of the Sun (Baraffe et al. 1998), however these models are 
not extended below O.GMq (My ~ 8.8). Baraffe et al. (1998) 
find that, for M < O.6M0, their isochrones with the three 
different mixing lengths converge, implying that the mix- 
ing length has little effect on the models below this mass. 
As such, we combine the a = 1.9 models for masses above 
~ O.6M0 with a = 1.0 models covering the very low mass 
stars. 

Lastly, we compare each group's method of deriving 
synthetic photometry for their models which is then used to 
convert effective temperatures to colours. The calculations 
by both the Geneva and Y^ groups come from versions of the 
BaSeL stellar spectral library (see eg. Lejeune and Schaerer 
2001). This library is based on model atmosphere spectra by 
Kurucz (1995), Fluks et al. (1994), and Bessell et al. (1989, 
1991). Padova calibrations are also based largely on the spec- 
tral library of Kurucz (1992) with extension to cool stars 
through the use of Fluks et al. (1994) empirical M giant spec- 
tra and the DUSTY99 (Chabrier et al. 2000) atmosphere 

^ http://pIeiadi.pd.astro.it/ lgirardi/isoc_photsys.html 
^ see reference for anonymous ftp instructions 
^ http://www.astro.yale.edu/demarque/yyiso.html 
^ http://webast.ast.obs-mip.fr/stellaj:/ 

http://www-laog.obs.ujf-grenoble.fr/activites/starevol/evol.html 



models of cool dwarf (typo M and later) synthetic spectra 
(sec §3 of Girardi et al. 2002 for a complete discussion). 
The Siess model uses Te// — BC relations from Schmidt- 
Kaler (1982) and Bessell (1991) and Teff — colour relation- 
ships from Fitzgerald (1970) and Arribas & Martinez- Rogers 
(1988, 1989). Lastly, the Baraffe isochrones use stellar spec- 
tra generated with the NextGen model atmosphere code (eg. 
Hauschildt, AUard, & Baron 1999) which is a predecessor of 
the aforementioned DUSTY99 models. 

In Fig. 4 we plot \og(L / Lsim) vs. log(re//) for each of 
the models used in this study (note that we have included 
the complete a = 1.0 as well as the a = 1.9 Baraffe models). 
All isochrones arc for solar metallicity and an age of 5 Gyr; 
the position of the Sun is indicated by the filled circle. The 
inset box shows a close-up of the region around the Sun. 
Comparing the models in the theoretical plane allows us to 
eliminate the uncertainty inherent in the transformation to 
the observational plane. In this figure wo see that four of 
the plotted isochrones, Padova (solid line), Baraffe, a = 1.9 
(dotted line), Y^ (dash-dot line), and Geneva (dash-dot-dot 
line) are all in good agreement for the upper MS and the po- 
sition of the Sun, but show considerable deviation beginning 
around O.ILq (Mv~7.5, O.65A/0, Teff ~ 4100K, and spec- 
tral type ~ K7). Neither Siess nor Baraffe (a = 1.0) match 
the Sun and, for a given luminosity, predict cooler effective 
temperatures than the other isochrones on the upper MS. 

Figure 5a displays the relationship between model mass 
and absolute magnitude in the V-band. In particular, the 
lower panel of Fig. 5a shows the difference in mass between 
each isochrone and the Padova set, which has been chosen 
as an arbitrary reference. We see that the mass differences 
are minimized at Mv~6.0, the magnitude at which we have 
determined the distance modulus of each cluster (Sec. 2). 
At fainter magnitudes the differences between the various 
models steadily increase but are generally less than O.IMq. 
Fainter than Mv~9, the masses of the models begin to 
deviate from the others. As we discuss how well the models 
fit the observations in the CMD, Fig. 5a will prove useful in 
estimating the mass at each value of the absolute magnitude. 

Figure 5b compares the difference in colour between 
each model and the Solar abundance Padova isochrone as a 
function of absolute V-band magnitude for {B — V)o (top 
panel), {V — I)o (middle), and {V — K)o (bottom). As 
with Fig. 5a, the colour differences are generally smaller 
at Mv^6-0 and increase at fainter magnitudes. Further- 
more, one would expect that if the models all used the same 
method for converting effective temperatures to colours, 
then the appearance of Fig. 5b would be solely a reflec- 
tion of the differences in the physics of the models (see 
Fig. 4). We see, however, that this is not the case; there 
arc no clear systematic trends in these differences from 
model to model (with the possible exception of the Geneva 
isochrones) . These virtually random differences indicate that 
the disagreement between models is a result of both the cho- 
sen physics and the colour transformations. 

4 RESULTS 

4.1 Color-Magnitude Diagrams 

As mentioned in §1, we axe interested in comparing the ob- 
servational CMDs with theoretical models for a variety of 
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colours. These comparisons are illustrated in Figs. 6 through 
22* for the 6 clusters in our sample. In each case, we show 
the [Mv,{B-V)o], [Mv,(y-/)o],and [Mv,{V-K)o] disr 
grams along with the isochrone that is closest to the clusters' 
mctallicity as listed in Table 1. Listed in each plot arc the 
cluster name and model name, along with two rmmbors that 
indicate which isochrone is used. The first rmmbor gives the 
log age of the model (e.g. 0945 \og{Age) = 9.45) and 
the second indicates the metal abundance (e.g. z019 =^ Z 
— 0.019). In the case of the Padova isochrones, Leo Girardi 
kindly provided interpolated models that match the clus- 
ter metal abundances. In addition, because the M37 data of 
Kalirai ct al. (2001) docs not possess I-band photometry, we 
only show data in the {B — V)o and {V — K)o colours. To 
simplify the discussion below, wo will refer to mismatches 
between the models and the observations as occuring in the 
magnitude direction with the understanding that there may 
be a component in the colour direction as well. 

Looking at the CMD and isochrone comparisons in Figs. 
6 through 22, a number of trends stand out. First, wc con- 
firm that all of the models fit the cluster MS at Mv=6.0 
in the {V — 7)o plane. This is to be expected given the 
manner in which wc determined the cluster distance moduli. 
Next, we examine the isochrone comparisons in the (_B — V^)o 
plane shown in Figs. 6 through 11. In the vast majority of 
cases, the model main sequences match the observations at 
Mv=6-0 adequately. There are two notable exceptions to 
this. First, the M37 photometry is significantly fainter than 
all of the isochrones considered herein by between ~0.5 and 
~0.9 mag in colour. As such, we will exclude it from further 
discussion within the context of the other {B — V)o CMDs. 
Second, the Baraffe models are consistently brighter than 
the observations at Mv=6.0, except in the case of M35. Mov- 
ing on to the other regions of the CMD in the {B—V)o plane, 
all of the models reproduce the unevolved MS in the range 
2.0^A4v^7.0 reasonably well. The primary exception to this 
is the fit of the Geneva models to the M35 photometry shown 
in the middle left panel of Fig. 6. In this case, the MS of the 
isochrone is clearly offset from the observed one for My ^5.0. 
Another exception, though not as obvious, is apparent in the 
comparison between the models and the M35 data in the 
upper right panel of Fig. 6. Again, the isochrone is brighter 
than the data in the range 2.0<C,My ^5.0. Thirdly, there is 
the comparison between NGC 2420 and the Siess models in 
Fig. 10 where the theoretical MS lies above the data. It is 
unclear what these mismatches arc due to. However, in the 
case of the Geneva models, it may be a result of the fact that 
the cluster (M35) metallicity is lower than those used in the 
isochrones. For the clusters whose photometry extends reli- 
ably below Mv~8.0 (NGC 2477 and M67, again excluding 
M37), we find that the Padova and Baraffe models diverge 
the most from the observed main sequences while the 
and Siess models deviate the least. Note that the Geneva 
model MS does not extend faint enough for this comparison 
to be made. 

Moving on to a consideration of the [My, {V — I)q] dia- 
grams in Figs. 12 through 16, we see a similar behaviour in 



* Only a sample of these CMDs (Figs. 6, 12, and 17) are included 
in the hard copy of this paper. The rest are only available in the 
electronic copy. 



the range 2.0^Afv<^7.0 as noted in the case of the {B — V)o 
CMDs. In short, all of the models perform reasonably well in 
matching the observed main sequences. The difficulties arise 
in the magnitude ranges fainter than Mv^S.O; the Padova 
models diverge the most from the cluster main sequences, 
while the Baraffe theoretical main sequences reproduce the 
observed ones more consistently than any of other sets. 

Lastly, we examine the [Mv, {V — K)o] diagrams which 
were constructed by matching optical photometry with 
2MASS ncar-IR photometry. Wo note that the level of agree- 
ment between the isochrones and the data at Mv=6-0 is not 
consistent from cluster to cluster and model to model. With 
the exception of M37 where the distance fits have been per- 
formed in the [My, {V — K)o] plane, in most cases, the model 
main sequences are fainter than the observed ones. However, 
the Padova MS matches that of M35 at My =6.0; similarly, 
the Geneva models match the MS of NGC 2477 at My =6.0. 
In terms of the level of overall agreement along the MS, we 
can say that, like the {B - V)o and {V - I)q CMDs, the 
{V — K)o models perform better in the range 2.0^Mv^7.0 
as compared to fainter magnitudes. However, with the ex- 
ception of M35, the models are generally fainter than the 
observational main sequences. 

To conclude this section, we note that none of the theo- 
retical models fit the unevolved MS (i.e. My ^4) adequately 
in all three colours. Either the shape of the isochrone does 
not match the observations or the location does not match 
or both. The small difference between the mctallicity of the 
isochrones and the observations is unlikely to be the cause 
of this mismatch. 



4.2 Two Color Diagrams 

One of the ways in which wc can eliminate the use of 
the distance modulus in the present analysis is to compare 
the optical/IR photometry of open clusters with theoretical 
isochrones in the two-colour diagram (TCD). In this plane, 
only the reddenings of the clusters are required allowing us 
to further investigate the robustness of the models. Figures 
23 through 28 show the [{V - K)o, {B - V)o\ TCDs while 
Figs. 29 through 33^ display the [{V - K)q, [v - I)o] TCDs 
for main sequence stars in our sample of open clusters. The 
arrow in each figure represents the slope of the reddening 
vector {E(B-V) = E{V-I) = 0.3) which is roughly parallel 
to the two-colour sequence suggesting that inferences made 
from these diagrams are largely insensitive to the adopted 
reddenings. 

Looking first at the [[V — K)o, [B — V)o] diagrams, ex- 
cept in a minority of cases, there is a zeropoint offset between 
the theoretical two-colour sequence and the observations; 
the models are too blue in {V — K)o at a given (B — V)o 
as compared with the observational data. In addition, the 
detailed shapes of the theoretical sequences can be dramat- 
ically different than the data especially in the case of the 
Baraffe models for {V - K)o^3.0. The M37 and M67 data, 
which extend to the reddest colours, nicely underscore these 
shape differences. Largely similar conclusions (i.e. models 

^ Again, only a sample of these figures are shown in the hard 
copy (Figs. 23 and 29) while the remainder are only available 
electronically. 
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arc offset to bluer colours and don't reproduce the shape 
of the observed TCDs) can be reached by examining the 
\{V - K)o, (V - I)o] diagrams. 

As a sanity check on the observations, Figs. 34 and 35 
present TCDs of each cluster compared with the sequence 
for M35 determined by eye using the data in the upper left 
panels. A reddening vector equal to E{B — V) = E{V — I) = 
0.3 is also shown in each figure. We see that in most cases, 
the positions and shapes of the two-colour sequences for each 
cluster arc in good accord with that of M35. In Fig. 34, the 
data for M37 are below the M35 line by ~0.1 mag; this 
may account for the fact that in the [My, {B — V)o] CMDs 
(Fig. 7), the M37 data is shifted away from the isochroncs 
by a significant amount as compared with the data from the 
other clusters. Figure 34 also shows that the M67 two-colour 
data appears to be slightly offset from the M35 sequence in 
the colour range 0.5^{B — y)o^l.O. Lastly, in Fig. 35, we 
see that the M67 sequence has a slightly different shape as 
compared with the M35 [{V - K)o,{V - I)o] line. Taken 
together, the TCDs in Fig. 34 and 35 provide reassurance 
that the observational data we have considered herein are 
internally consistent. 



5 CONCLUSION 

We have presented combined optical/near-IR photometry 
for six open clusters. The optical data have been obtained 
as part of the WIYN Open Cluster Study and taken from 
the literature, while the near-IR observations are from the 
Second Incremental Data Release of the Two Micron All Sky 
Survey Point Source Catalog. The open clusters span an age 
range from 150 Myr to 4 Gyr and have metal abundances 
from [Fe/H] = -0.27 to -1-0.09 dex. We have utihzed these 
data to test the robustness of theoretical main sequences 
constructed by several groups as denoted by the following 
designations - Padova (Girardi et al. 2002), Baraffe (Baraffe 
et al. 1998), (Yi et al. 2001), Geneva (Lejeune & Schaerer 
2001), and Siess (Siess et al. 2000). The comparisons of the 
models with the observations have been performed in the 
[Mv, {B - V)o], [Mv, {V - /)(,], and [My, (V - K)o] colour- 
magnitude diagrams as well £is the distance-independent 
[{V - K)o, {B - V)o] and [{V - K)o, (V - I)o] two-colour di- 
agrams. The reddenings have been talcen from Twarog et al. 
(1997) and the distance moduli were determined by adopting 
these reddenings and then matching each isochrone to the 
observed main sequences at Mv=6-0 in the [My, {V — I)o] 
plane {[My, {V - K)o] for M37). 

We conclude that none of the theoretical models repro- 
duce the observational data in a consistent manner over the 
magnitude and colour range of the unevolved main sequence. 
In particular, there are significant zeropoint and shape dif- 
ferences between the models and the observations. In the 
colour-magnitude diagrams, most of the models performed 
reasonably well for stars with MassJ>O.75M0, but did not 
do so for stars with lower masses. In the two-colour dia- 
grams, the theoretical sequences are too blue in (V — K)o at 
a given {B — V)o and (V — /)o as compared with the observa- 
tional data. The shapes of the theoretical sequences can also 
be dramatically different as compared with the data. From 
Figs. 4 and 5 we can say with a fair amount of certainty that 
the crux of the problem lies in the precise mismatch between 



theoretical and observational colour-temperature relations. 
Depending on how these relations are derived, one of the 
problems could be missing sources of opacity in the atmo- 
spheres of the model stars. These results underscore the im- 
portance of pursuing the study of stellar structure and stellar 
modelling with even greater intensity. 
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Figure 1. Colour-magnitude diagrams in the [Mv,{B — l/)o] 
plane for the six clusters in our dataset. The distance moduli and 
reddenings are given in Table 1. All plots have the same axis scale 
and the clusters axe presented in order from youngest (top left) 
to oldest (lower right). 



Figure 2. Same as Fig. 1 except that the CMDs are in the 
{Mv,{V — /)o] plane. We note that M 37 is not included be- 
cause Kalirai et al. (2001) do not include I-band photometry in 
their dataset. 



Figure 3. Same as Fig. 1 except that the CMDs arc in the 
[My, {V — K)q\ plane. It is evident from the plots that the CMDs 
are limited by the depth of the K-band data. 

Figure 4. We plot log(L/L0) vs. \og,{T^ff) for each of the models 
used in this study. All isochrones are for solar metallicity at an 
age of 5 Gyr. The position of the Sun is marked by the filled 
circle. The inset shows a close-up of the region around the Sun's 
location. 

Figure 5. (a) The upper panel shows the relationship between 
absolute magnitude in the V-band and stellar mass in units of 
solar masses for each of the models considered herein. In the lower 
panel, we have arbitrarily selected the Padova isochrones as the 
reference and plot the mass difference between eaich set of models 
and this reference. 



Figure 5. (b) In these three ijaiicls wc plot the difference in colour 
between each model and the Solar abundance Padova isochrones 
as a function of absolute V-band magnitude for (B — V)o (top 
panel), {V — I)o (middle), and {V — K)o (bottom). 



Figure 6. Comparisons between the [Mv,{B — V)o] CMD of 
M35 and theoretical models from the Padova, Y^, Geneva and 
Siess groups. Note that the cluster distance modulus has been 
adjusted so that the models and observations agree at Mv^=6.0 
in the VI CMDs. The Log of the plotted ages are given along with 
the metal abundance (Z) of the models. 



Figure 7. Same as Fig. 6 except that the photometry of M37 is 
shown. 



Figure 8. Same as Fig. 6 except that the photometry of NGC 
1817 is shown. 



Figure 9. Same as Fig. 6 except that the photometry of NGC 
2477 is shown. 



Figure 10. Same as Fig. 6 except that the photometry of NGC 
2420 is shown. 

Figure 11. Same as Fig. 6 except that the photometry of M67 
is shown. 



Figure 12. Comparisons between the [My, {V — I)o] CMD of 
M35 and theoretical models from the Padova, , Geneva, Siess, 
and Bajraffe groups. The Log of the plotted ages are given along 
with the metal abundance (Z) of the models. 
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Figure 13. Same as Fig. 12 except that the photometry of NGC 
1817 is shown. 



Figure 14. Same as Fig. 12 except that the photometry of NGC 
2477 is shown. 



Figure 15. Same as Fig. 12 except that the photometry of NGC 
2420 is shown. 

Figure 16. Same as Fig. 12 except that the photometry of M67 
is shown. 

Figure 17. Comparisons between the [My, (V — K)q] CMD of 
M35 and theoretical models from the Padova, Y^, Geneva, Siess, 
and Baraffe groups. The Log of the plotted ages are given along 
with the metal abundance (Z) of the models. 

Figure 18. Same as Fig. 17 except that the photometry of M37 
is shown. 

Figure 19. Same as Fig. 17 except that the photometry of NGC 
1817 is shown. 



Figure 20. Same as Fig. 17 except that the photometry of NGC 
2477 is shown. 



Figure 21. Same as Fig. 17 except that the photometry of NGC 
2420 is shown. 

Figure 22. Same as Fig. 17 except that the photometry of M67 
is shown. 

Figure 23. Two colour diagrams (TCD) using {B — V)o and {V — 

K)o data for main sequence stars in M35. The solid lines represent 

the models indicated in each panel. These plots show that even 
when we remove the uncertainties inherent in the distance moduli, 
the theoretical isochroncs still do not reijroduce the shape of the 
observed sequence in the TCD. We have included a reddening 
vector with length {B - V)o = 0.3 



Figure 24. Same as Fig. 23 except that we now compare the 
observational data for M37 with theoretical models. 
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Figure 25. Same as Fig. 23 except that we now compajre the 
observational data for NGC 1817 with theoretical models. 



Figure 26. Same as Fig. 23 except that we now compare the 
observational data for NGC 2477 with theoretical models. 



Figure 27. Same as Fig. 23 except that we now coiripare the 
observational data for NGC 2420 with theoretical models. 



Figure 28. Same as Fig. 23 except that we now compare the 
observational data for M67 with theoretical models. 



Figure 29. Same as Fig. 23 except that two colour diagrams 
using (V — /)() and {V — K)o data are plotted. The reddening 
vector shown has length {V — I)o = 0.3. 



Figure 30. Same as Fig. 29 except that we now compare the 
observational data for NGC 1817 with theoretical models. 



Figure 31. Same as Fig. 29 except that we now compare the 
observational data for NGC 2477 with theoretical models. 



Figure 32. Same as Fig. 29 except that we now compare the 
observational data for NGC 2420 with theoretical models. 



Figure 33. Same as Fig. 29 except that we now compare the 
observational data for M67 with theoretical models. 



Figure 34. Two colour diagrams using {B — V)o and {V — K)(i 
data for all of the clusters in the present study compared with 
the observed sequence for M35 (solid line). We have included a 
reddening vector with length {B — V)o = 0.3 



Figure 35. Two colour diagrams using (V — /)<) and {V — K)o 
data for all of the clusters in the present study compared with 
the observed sequence for M35 (solid line). We have included a 
reddening vector with length {V — I)o = 0.3 



Figure 36. The [My, {V - K)o] CMDs for M35 and M67 with 
the solid lines representing the fiducial sequences drawn by eye 
through the highest density of points. The right panel shows the 
fiducials for both clusters plotted in the same plane along with 
field stars (open circles) from Percival et al. (2003) with close-to- 
solar abundances and Hipparcos parallaxes. 
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Table 1. Ojjeii Chistor Iiiforiiiatioii 



Name 


Available Photometry 


Log Age 


E{B - V) 


[Fe/H] 


M 35 (NGC 2168) 


UBVRIJHKs 


8.17 


0.19 


-0.160 


M 37 (NGC 2099) 


...BV JHKs 


8.73 


0.27 


0.089 


NGC 1817 


...BVRIJHKs 


8.80 


0.26 


-0.268 


NGC 2477 


UBV...IJHKs 


9.04 


0.23 


0.019 


NGC 2420 


...BVRIJHKs 


9.24 


0.05 


-0.266 


M 67 (NGC 2682) 


UBVRIJHKs 


9.60 


0.04 


0.000 



Table 2. Open Cluster Distance Moduli. 



Cluster 


Padova 


Baraffe 


Geneva 


Y2 


Siess Twarog et al. 


M 35 (NGC 2168) 


10.16 


10.41 


9.81 


9.91 


9.96 


10.30 


M 37 (NGC 2099) 


11.55 


11.40 


11.50 


11.35 


11.75 


11.55 


NGC 1817 


12.10 


12.30 


11.90 


11.85 


12.00 


12.15 


NGC 2477 


11.55 


11.60 


11.30 


11.15 


11.45 


11.55 


NGC 2420 


12.12 


12.45 


11.95 


11.90 


12.07 


12.10 


M 67 (NGC 2682) 


9.80 


9.80 


9.60 


9.45 


9.65 


9.80 
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Table 3. Theoretical Model Input Parameters 



Padova 



Baraffe 



Geneva 



Y2 



Siess 



Opacity 

Low temperature 

opacity 

Equation of state 



Core overshoot 

Mixing length, a 
He abundance 
He enrichment, 

Synthetic photometry 



OPAL (1993)" 
AF94 



T > 10'': Kippenhahn"* 
T < 10'^: MHD" 



M ^ 1.5M0 
1.68 
Yp = 0.23 
2.25 



ATLASQf 
DUSTY99'' 
Fluks et al. (1994) 



OPAL (1996)'' OPAL (?)= 

AF94 AF94, 

Kurucz (1991) 

SCVH-'' Maeder & 

Meynet (1989) 



M ^ 1.5M0 
1.9 1.6 
Y,„iar = 0.282 Yp = 0.24 



n/a 



NextGen* 



2.5 for Z > 0.02 
3 for Z ^ 0.02 

BaSeL-2.2J' 



OPAL (1996)'' OPAL (1996)*' 
AF94 AF94 



OPAL (1996)'' based on 

Pols et al. (1995) 



o-ge <, 2Gyr (all others = 0) 

1.7431 1.605 



2.0 



2.0 



BaSeL-2.2-' Siess et al. (1997) 



"Iglcsias, & Rogers (1993) 
''Iglesias, & Rogers (1996) 

^Geneva isochrones were published over the course of several years 

and as such utilize OPAL opacities from different years. See Leje- 

une & Schaerer (2001) for more information. 

'^Kippenhahan, Thomas &; Weigert (1965) 

''Mihalas et al. (1990) 

^Saumon, Chabricr, & VanHorn (1995) 

sCastelli, Gratton, & Kurucz (1997) 

'"Chabricr ct al. (2000) 

'Hauschildt, AUard, & Baron (1999) 

^Westera, Lejeune, & Buser (1999) 
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Table 4. Cluster Fiduclals locations of the field stars and the MS locations of M35 and 

' M67. 

M35 M67 



V 


V - K 


V 


V - K 


9.0 


0.32 


11.0 


2.57 


9.5 


0.34 


11.5 


2.44 


10.0 


0.36 


12.0 


2.32 


10.5 


0.40 


12.5 


2.20 


11.0 


0.44 


12.9 


2.09 


11.5 


0.51 


13.0 


2.01 


12.0 


0.66 


12.8 


1.85 


12.5 


0.89 


12.7 


1.74 


13.0 


1.17 


12.65 


1.58 


13.5 


1.45 


12.7 


1.40 


14.0 


1.66 


12.8 


1.33 


14.5 


1.83 


13.1 


1.42 


15.0 


2.01 


13.2 


1.41 


15.5 


2.19 


13.4 


1.38 


16.0 


2.39 


13.5 


1.375 


16.5 


2.64 


13.6 


1.374 


17.0 


2.96 


13.7 


1.385 


17.5 


3.29 


13.8 


1.40 


18.0 


3.59 


14.0 


1.43 


18.5 


3.89 


14.2 


1.48 


19.0 


4.19 


14.4 


1.55 


19.5 


4.44 


14.6 


1.62 


20.0 


4.69 


14.8 


1.69 






15.0 


1.77 






15.5 


1.99 






16.0 


2.26 






16.5 


2.54 






17.0 


2.84 






17.5 


3.17 






18.0 


3.50 






18.5 


3.80 






19.0 


4.07 



APPENDIX A: CLUSTER V -K FIDUCIAL 
SEQUENCES 

Published {V — K)o fiducial sequences are rare for open clus- 
ters; as such, to help fill an important need, we have derived 
such sequences for the youngest (M35) and oldest (M67) 
clusters in our sample. The left and centre panels of Fig. 36 
show our adopted fiducial sequences for M35 and M67, re- 
spectively. These have been constructed by-eye through the 
highest density of points iu the CMD and arc given in Table 
4. The cluster data and the fiducials are plotted after ac- 
counting for the effects of reddening (Table 1) and distance 
using the median of the values in Table 2. The right panel in 
Fig. 36 shows the two fiducials plotted together to reinforce 
the point that the observations of these clusters agree well 
with each other even to the limits of the M67 photometry 
at Mv^Q. 

To investigate how well solar-abundance field MS stars 
compare with these two clusters, we have combined optical 
photometry of field stars with [Fe/H] values within ±0.1 
of the Sun and Hipparcos parallaxes from Table 1 of Perci- 
val, Salaris, & Kilkenny (2003) with K-band data for these 
same stars from the Second Incremental Data Release of the 
2MASS Point Source Catalog. The open circles in the right 
panel of Fig. 36 represent three stars that meet these crite- 
ria. It is evident that there is good agreement between the 
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